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ABSTRACT 


Electron density deviations from a basic variation with the 
solar zenith angle are investigated. A model study has been 
conducted in which the effects of changes in neutral and relative 
densities of atomic and molecular oxygen on calculated electron 
densities have been compared with incoherent scatter measurements 
in the height range 100- 117 km at Arecibo, Puerto Rico. The 
feasibility of determining tides in the neutral atmosphere from electron 
density profiles was studied. It was determined that variations in 
phase between the density and temperature variation and the 
comparable magnitudes of their components makes it appear 
improbable that the useful information on tidal modes can be obtained 
in this way. 


CHAPTER I 


GENERAL STATEMENT OF THE PROBLEM 
1.1 E-Region Morphology 

Appleton (1953) discussed the differences in time between 
the maximum of ionization and the maximum of production in terms 
of the "sluggishness" of the E-region. Researchers (Mitra, 1958; 
Appleton and Lyon, 1961; and Palluconi, 1963) attempted to 
determine the recombination rate, a, of the E-regio' by observing 
asymmetries about noon in foE using the charge continuity equation. 
The lag or sometimes lead value t was employed. While many 
attempts were made to calculate recombination rates in this way it 
became apparent with modern laboratory measurements (Biondi, 
1964) that the time constants were too small to account for the 
effect and that the explanation of the phase differences must be else- 
where. 

Monro, Nisbet and Stick (1976) in the first part of the present 
study demonstrated that local noon asymmetries in the E- region 
electron density were caused by tides in temperature, density, and 
winds as measured by incoherent scatter. They demonstrated that 
the effect of the semidiurnal tide on the atomic and molecular oxygen 
densities may also be an important factor in explaining E-region 
asymmetries. 

These studies can only be done at a very few locations none 
of which are in the southern hemisphere. The global tidal variations 
are extremely important in determining the mode structure and 
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wider geographical coverage is necessary before the tidal morphology 
can be resolved. The tidal variations in the lower thermosphere have 
been shown to be very variable (Wand and Perkins, 1970; Salah and 
Wand, 1974; and Fellous et al. , 1975) and thus a long series of 
observations is desirable. It, therefore seemed important to 
reexamine the studies that had been made of E- region changes to 
determine to what extent the data can now be used to study the tides. 

A complete analysis of the electron density at fixed heights is 
necessary to determine morphology of the tidal system. Some 
indication may, however, be determined from the one parameter 
that has been extensively studied, t, the difference in time between 
the minimum solar zenith angle and the maximum electron density. 
Palluconi (1963) has shown that t varies with latitude such that it is 
generally negative for equatorial stations becoming positive at 20° 
and reaching a peak at 40^. The value of t is generally symmetric 
about the equator during equinox with a shift in the minimum from 
the equator toward the winter pole during solstic?. The general 
seasonal response for high latitude stations is a maximum value of 
t in the summer with a minimum in the winter and generally a 
n^gati^e value throughout the year for equatorial stations. Since the 
number of southern hemisphere stations is minimal it is not possible 
to conclude with great certainty the tidal effects simultaneously in 
both hemispheres from previous studies. Data is available in the 
World Data Center for several years that could, however, be 
analyzed. 


1 . 2 Solar Fluxes 


The solar wavelengths of importance in photoionizing the 
E-region have been reviewed by Rishbeth and Garriott (1969) and 
Banks and Kockarts (1973). X-rays from 8 to 140 X are important 
in ionizing the three major constituents, molecular nitrogen and 
atomic and molecular oxygen in the lower E-region, particularly in 
times of enhanced solar activity. The EUV flux from 140 to 78t'> X 
is important in ionizing molecular nitrogen and atomic oxygen in the 
upper E-region. Wavelengths from 786 X to 1027 X are important 
with molecular nitrogen being ionized up to 796 X, atomic oxygen 
up to 911 X and molecular oxygen up to 1027 X including Lyman P. 
Nitric oxide is a very minor constituent in the upper atmosphere 
however as was shown by Nicolet and Aitkin (1959) it is of importance 
in the D-region and lower E-region because of its aoility to be 
photoionized up to 1340 X which includes H Lyman a. 

Donnelly and Pope (1973) compiled the results of different 
solar flux measurements and normalized them to a 10.7 cm flux of 
150 x 10"^ Wm*^ Hz" * at 1 AU. Their solar flux model includes 
wavelengths of 1 to 3000 X. As 3hown by Hall and Hinteregger (1970) 
care must be ujed when applying a flux model for a different 10.7 cm 
flux value then for the measurement conditions. Table 1 shows the 


solar flux values of Hinteregger et al., (1965) for low solar activity, 
^HHS’ ant * °* Donnell y anc * P°P e (1973) for medium solar activity, 


Variations in EUV spectral lines over one solar rotation 


period of 27 days have been reported by Hall and Hinteregger (1970) 
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Absorption Cross Sections, Photoionization Yields and Assumed Photon Fluxes 



Wavelength 

Absorption Cross Sections 

Photoionization Yields 

o 

Photon Fluxes $/10 

, - 2 

Ph-cm -sec 

No. 

\ X 

/in' 18 2 

or / 1 0 cm 










o 

O- 

N, 

o 

o. 

N 

$ 

<t> 





2 

2 


2 

2 

HHS 

DP 

DP 

1 

15-10 

. 135 

.27 

.15 

16 7 

16.7 

18.9 

.0005 

. 0092 

16-10 

2 

22.8-15 

.35 

.7 

.36 

12.4 

12.4 

13.8 

.003 

. 015 

23.3-16 

3 

31-22.8 

.045 

.09 

1.0 

10. 1 

10. 1 

12.8 

.005 

.015 

31-23.3 

4 

41-31 

. 1 

.2 

.07 

8.4 

8.4 

9.4 

.02 

.036 

40-31 

5 

50-41 

. 17 

. 34 

. 14 

6.7 

6.7 

7.6 

. 03 

. 040 

50-40 

6 

60-50 

.27 

. 54 

.22 

6.7 

6.7 

7.6 

.068 

.084 

60-50 

7 

70-60 

.35 

.7 

.32 

4.8 

4.8 

5.4 

. 08 

.099 

70-60 

8 

80-70 

.45 

.9 

.42 

4.8 

4.8 

5.4 

. 08 

.094 

80-70 

9 

90-80 

. 5 

1.0 

.4 

4.2 

4.2 

4. 8 

. 1 

. 13 

90-80 

10 

103-90 

.9 

1.8 

.7 

3.7 

3.7 

4.2 

. 13 

. 106 

100-90 

11 

139-103 

1.1 

2.2 

.9 

2.9 

2.9 

3.3 

.304 

. 153 

140-100 

12 

165-138 

2. 1 

4.2 

1.4 

2.3 

2.3 

2.6 

.6 

.32 

170- 140 

13 

180-165 

2.9 

5.7 

1 9 

2.0 

2.0 

2.3 

.5 

1.47 

180-17U 

14 

190-180 

3.9 

6.8 

2.2 

1.9 

1.9 

2. 1 

2. 2 

1.15 

190-180 

15 

205-190 

4.0 

8.0 

2.4 

1.8 

1.8 

2. 0 

1.5 

1.4 

210-190 

16 

220-205 

4.7 

9.4 

2.3 

1.6 

1.6 

1.6 

.55 

.31 

220-210 

17 

240- 22C 

5.6 

11.2 

3. 1 

1.2 

1.2 

1.2 

.85 

1.34 

240-220 

18 

260-240 

6.7 

13.4 

3.6 

1.2 

1.2 

1.2 

1.3 

1.23 

260-240 

19 

256.3-257 

7.2 

14.4 

3.8 

1.2 

1.2 

1.2 

.40 

. 59 

256. 3 

20 

280-260 

8.0 

16.0 

4.2 

1.0 

1.0 

1.0 

1.4 

1.08 

280-260 

21 

810-796 

3. 3 

40.0 

6.7 

1.0 

2.4 

0 . 0 

.2 

.25 

810-796 

22 

840-810 

3.2 

26.0 

3.3 

1.0 

.38 

0 . 0 

.78 

1.00 

840-810 

23 

832-835 

3.2 

13.5 

11.1 

1.0 

.32 

0.0 

.52 

. 53 

834. 2 

24 

860-840 

3. 1 

11.7 

5. 0 

1.0 

.31 

0.0 

i . 5 

1.09 

860-840 

25 

890-860 

2.9 

8.9 

5.0 

1.0 

.63 

0 . 0 

2.7 

2. 78 

890-860 

26 

911-890 

2.7 

9.3 

4.3 

1.0 

.75 

0 . 0 

3.33 

3.23 

911-890 

27 

920-91 1 

0.0 

7.4 

3.8 

0.0 

.65 

0.0 

1.25 

2. 30 

920-911 



Table 1: Absorption Cross Sections, Photoionization Yields and Assumed Photon Fluxes (continued) 
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compared with the 10.7 cm flux which ranged from 100-200 x 10 

Wm’^ Hz'*. The ratio of the fluxes 284 R to 304 R was shown to 

increase by 200 per cent as the 10.7 cm flux increased from 100 to 

200 x 10 ^ Wm ^ Hz *. The combined fluxes of 1025 A°, 977 R, 

584 R and 304 R increased by 70 per cent as the 10.7 cm flux 

- 22 - 2 - 1 

increased from 60 to 160 x 10 Wm Hz 

A discrepancy has been shown to exist between the results of 
ion densities calculated using current neutral atmospheric models 
and photochemistry with the fluxes of Hinteregger (1970) and 
observations of ionosphere parameters. Swartz and Nisbet (1973) 
comparing profiles of energy input and loss with similar profiles 
of ion production and loss concluded that the EUV fluxes needed to be 
multiplied by a factor of two. Mitra et al. , (1974) found the fluxes 
of soft x-rays to be incompatible with the observations of electron 
density and the ratio of M 1 [ 0 + ] in the altitude range of 90 to 
130 km. It was shown that an increase by a factor of three for the 
soft x-rays is needed to reconcile calculations and observations. 
Heroux et al. , (1974) modeled the region of 110 to 300 km and found 
an electron density that was generally low by 30 per cent compared 
with an electron density profile derived from an ionogram. They 
concluded that the error was probably in the density measurements. 
Nisbet (1975) showed that a general doubling of the flux would increase 
the electron densities to match the experimental values. Heroux et 
al. , (1975) replied that while a doubling of the flux would help at the 
lower altitudes the increase at the upper altitudes would be too great. 
Donahue et al., (1973) and Monro, Nisbet and Stick (1976) have 
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shown that the atomic oxygen is very variable and can differ greatly 
from models such as CIRA 1972. It is difficult to obtain absolute 
values for solar fluxes corresponding to given solar conditions. The 
discrepancies appear to be greater under high solar conditions than 
for low activity (Miller and Nisbet, 1974). All of these factors 
contribute to making such comparisons difficult. 

Ionization cross sections were reviewed by Ohshio et al. , 
(1966). The absorption cross section along with the photoionization 
yield for the major E- region neutrals, N^, and O corresponding 
to different regions of the EUV and x-ray spectrum are listed in 
Table 1. Swider (1969) described the production of secondaries by 
x-rays as a function of the total production, the density of the 
constituents and the wavelength intervals. The total ion- electron 
pair production rate q was shown to be distributed among the 
different ions according to the ratio 0.62 : 0.17 : 0.14 : 0.07 : 1.00 
for q(N^) : q(0*) : q(N + ) : q(0 + ) : q total . 

1.3 E-Region Tides 

Absorption of solar energy by ozone in the strato-mesosphere 
causes waves to propagate into the thermosphere. Chapman and 
Lindzen (1970) modeled the atmosphere and derived variations in 
temperature, wind and density in terms of altitude, latitude and 
local time. Richmond (1971) included dissipative effects which 
resulted in a dampening of the diurnal (1,1) mode above 105 km 
while the semidiurnal modes were free to propagate through the 
region. Bernard and Spizzichino (1971) using incoherent scatter 
techniques were able to show that the semidiurnal tide dominated 
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both wind and temperature oscillations above 100 km. At altitudes 
below 100 km the diurnal tide, the semidiurnal tide and gravity 
waves have comparable amounts of energy. Salah (1974) observed 
a dominant semidiurnal tide in the lower thermosphere in both 
temperature and meridional winds corresponding to the (2,4) model. 
Lindzen and Hong (1974) developed a program which showed a great 
amount of mode coupling between the main semidiurnal mode (2, 2) 
and higher modes which results in higher modes of the semidiurnal 
tide dominating above 100 km and a reduction in the amplitude of the 
(2, 2) mode. 

Salah et al., (1975) discussed results from three incoherent 
scatter stations for the lower E- region. A good agreement in both 
magnitude and phase of the average semidiurnal tides was found for 
Millstone Hill and St. Santin. The tide was identified as a (2,4) mode. 
The oscillations at Arecibo were dominated by the semidiurnal tide 
however great variability in phase was observed. The implications 
of simultaneous experiments at the three stations are that the 
E-region is controlled by local effects. 

Seasonal variations of the tides were observed by Salah and 
Wand (1974) for Millstone Hill and St. Santin revealing a smaller 
amplitude in the semidiurnal tide during the winter than during the 
summer or equinox due to a larger phase variability in the winter. 

Salah (1974) using incoherent scatter data showed the neutral 
density variations to be semidiurnal and in phase with the temperature 
at the lower E- region altitudes. The variations for the neutral 
density were of the order of a factor of two while temperature varied 
by 1 1 per cent. 
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Mayr and Vollan*! (1972) used a neutral atmospheric model to 
predict that the atomic oxygen density would vary diurnally with a 
maximum shifted toward early morning by the zonal wind and toward 
later afternoon by vertical diffusion in the daytime. The net result 
is a maximum at 1030 local solar time. Alcayde et al., (1972) 
using measurements of satellite drag and incoherent scatter show 
atomic oxygen variations at 120 km, extrapolated down from higher 
altitudes, to have a dominant diurnal mode with a superimposed 
semidiurnal mode. There are maximums at 900 and 1500 hours with 
a minimum at 1200. The ratio of maximum to minimum is 
approximately 1.5. Danilov (1972) using rocket mass spectrometer 
data derived a model of atomic oxygen starting at 130 km for low, 
moderate and high solar activity. For all three solar conditions 
the model shows a diurnal variation in atomic oxygen with a maximum 
at 1300 hours. During moderate activity a ratio of two is reported 
for maximum to minimum density at 130 km. Offerman (1974) 
using rocket measurements finds a maximum in molecular 
nitrogen to occur at 1400 hours at 150 km. Molecular oxygen was 
found to have a diurnal variation with a maximum at 900 to 1000 hours 
with a 20 per cent variation. 

1.4 Ion Chemistry 

The major neutral constituents photoionized in the E- region 
are N^, and O. Nitric oxide ionized by H Lyman a is also of 
importance. The photochemistry of the E- region has been reviewed 
by Biondi (1969), Bortner et al. , (1972) and Ferguson (1974). The 
rates used in this work have been determined from a review of current 
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literature. They are presented in Table 2. 

1.5 Sporadic E and Tides 

Whitehead (1961) explained sporadic-E formation as being due 
to vertical wind shears in the horizontal neutral winds. Whitehead 
(1967) elaborated cn his original explanation by demonstrating that 
the existence of a second ion with a smaller recombination rate, such 
as a metallic ion, could explain the general shape of the sporadic-E 
layer. The theory showed that the vertical wind shear of the zonal 
wind was more important than that of the meridional in the lower 
E-region for forming a sporadic-E layer. 

Wright et al. , (1967) observed neutral winds from visible 
trails of gun launched vehicles in conjunction with radio soundings 
of the sporadic-E layer. Their results show that sporadic-E occurs 
where the vertical ion movements converge. The wind profile was 
observed to descend at the rate of two vertical wavelengths in 24 
hours . 

MacDougall (1974) used the theory of wind shear to determine 
characteristics of the neutral zonal winds in the lower E- region 
from ionogram measurements of the sporadic-E layer altitude and 
it's descent with time. He assumed that the sporadic-E layer formed 
where the neutral wind above had a westward direction and the wind 
below an eastward direction for layers existing below 115 km. Above 
about 120 km the vertical shear in the north- south neutral wind 
determined the convergence of ionization. The ionogram data on 
the sporadic-E layer height evolution revealed three general patterns; 



Table 2: Chemical Reaction Rates (continued) 
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Table 2: Chemical Reaction Rates 
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a semidiurnal, a diurnal, and a steady wind component. The semi- 
diurnal wind component predominated in the middle latitudes while 
being spai se in the equatorial zone and at high latitudes. The results 
indicated a mixture of several semidiurnal modes. The semidiurnal 
components were found to be in good agreement with that of Bernard 
and Spizzichino (1971). The diurnal component was found mainly for 
equatorial stations and weak in the middle latitudes. The steady 
wind component was observed to be generally eastward in a band of 
± 10° of the equator and generally westward in the region of 10° - 30° 
from the eauator. 

Harper et al. , (1975) found the spcradic-E layer to move in 
accord with the convergence node in the vertical ion motion. Above 
124 km the derived wind in the north-south direction had a node 
near the sporadic-E layer with a poleward wind above and an equatorial 
wind below. 

Nighttime E-region ionization sources have been reported 
by Young et al. , (1971), Fujitaka et al. , (197i), Johnson (1972) and 
Strobel (1974) to be primarily H Lyman 8 interplanetary emission 
of He I 504 R and He II 304 R, starlight from 1026 to 910 R and a 
miner contribution from H Lyman a. Fujitaka et al. (1971) observed 
peak electron density layers in conjunction with rocket measurements 
of vertical ion velocity. The layers of peak density were found to 
occur where the ior.s converged. The points of convergence were 
found to descend througnout the night. Geller et al., (1975) using 
nighttime rocket measurements of electron der.s.ty in the continuity 
equation were able to reconstruct the vertical ion drifts. The drift 
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velocity was a function of the shape of the enhanced ionization layer. 
The ion drift so determined indicated a tidal wind of the (2,4)mode. 

1.6 Specific Statements of the Problem 

1. To investigate the effects of variations in neutral density, 
vertical ion velocity and temperature on the E- region electron density 
profiles and to determine the importance of each in making midlatitude 
electron density models. 

2. To study effects of tides in the E-region on the electron 
densities and to examine the feasibility of determining tidal 
information from electron density profiles. 

3. To compare measured solar fluxes and electron ion 
densities with that of a model and to examine causes of any 
discrepancies found. 


CHAPTER II 


E-REGION MODEL 

2. 1 Neutral Model Atmosphere 

The E-region is considered to be that part of the ionosphere 
immediately above the turbopause and tnerefore its' major neutral 
constituents may be assumed to follow a molecular diffusion profile. 
As such the altitude profiles of the major constituents O, O^ and N^ 
are determined by the individual scale heights 


H. 

1 


kT 
m. g 


( 2 . 1 ) 


where 

k is Boltzmann's constant 

T is the neutral temperature 

g is the acceleration due to gravity 

th 

rrr is the molecular weight of tne i constituent 
having the height relation for element n. 


n.(z) = n. T /T(z) exp 

i i o r 

o 

where n. and T are the number density and lemperature respectively 
l o ° 

of constituent i at altitude z q and T(z) is the temperature at altitude z. 
For the initial model it was assumed that the relative densities of O, 
and N^ at each height was as given in CIRA (1972). The density 


z 

i 


dz 

H.(z) 


( 2 . 2 ) 
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of the minor constituent nitric oxide was taken to be the model values 
reported by Meira (1971). These were in Bom. cases modified as 
will be described later. The total amount of any one constituent 
above the considered altitude can be approximated by 


00 

J n. dz = n. H. 

1 l o 1 

o 


(2.3) 


Values of ion collision frequency and ion temperature were 
obtained by analyzing the incoherent scatter data from Arecibo, 
Puerto Rico for the six altitudes 100, 104, 107, 110, 114 and 117 km. 
The total neutral density was derived irom the ion collision frequency 
by the relationship of Banks (1966). Because the accuracy of the 
temperature measurements is greater than the accuracy of the 
collision frequency measurements at the upper and lower ends of 
the altitude range of interest, an averaging procedure based on 
diffusion equilibrium was used. Diffusion equilibrium can be 
assumed since the vertical accelerations induced by the vertical 
motiors, including the tides, are much smaller than that due to 
gravity. It is thus much better to use all the collision frequencies 
measured to derive a best estimate of the collision frequency of each 
measurement and the corresponding temperature. This average 
collision frequency can then be used to derive the neutral densities 
at the six altitudes using scale heights calculated from the 
temperature measurements. In this work the collision frequencies 
were normalized to 110 km according to the formula 
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(2.4) 

where is the collision frequency at 110 km based on the j 1 * 1 

altitude and is the measured collision frequency at the j tl1 altitude. 
The average v jjq can then be calaculated as 



(2.5) 


where is the error in collision frequency at the height. 

The measured ion temperatures were assumed to be equal 
to the neutral temperatures as shown by Nisbet (1967) for E-region 
altitudes. As the accuracy of the temperature measurements was 
considerably better than that of the ion collision frequencies, the 
reduced temperature measurements were adopted for each height. 


2.2 Solar Effects and Photoionization 

The equation for calculating the production of ionization by 
solar photons can be derived using the well known equation (2.6), 
see for example Rishbeth and Garriott (1969), 


Q(\, z) 


£ $ exp 

. i CD r 


-sec y 


£ a. . H.(z) £ V, n (z) 

j ‘J J J j l J ‘J J 


( 2 . 6 ) 
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where the subscripts i and j refer to the wavelength group and the 
individual constituents respectively and 




is the photon flux outside the earth's atmosphere at one 
astronomical unit, 
is the solar zenith angle. 

is the total absorption cross section for constituent j 
and wavelength i. 

is the scale height of the j 1 * 1 constituent. 

is the photoionization yield for constituent j and 

wavelength group i. 


n.(z) 


is the density of the constituent. 



in Table 1. The individual neutral densities and scale heights were 
computed as explained in section 2. 1. The photon flux for each 
wavelength gioup was taken to be the Hinteregger et al., (1965) value 
repeated on Table 1. The angle X was calculated as a function of 
time of day, latitude and solar declination. 


X = Arccosine sin (g) sin (6) + cos (g) sin (6) cos | ~ * ^ j 


(2.7) 


where 

g is the geographical latitude. 

HL is the local solar time in decimal hours. 
6 is the solar declination angle. 
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In this study it has been assumed that 


6 


- .40915 


cos 


2 tr (D + 8) 
365. 25 


radians 


( 2 . 8 ) 


where D is the day number. 

The effect of secondaries being formed by x-rays was modeled 
according to the method described by Swider (1969). First the 
total production was determined as stated above accounting for 
increased production through the ionization yield. By energetic 
photoelectrons the ion production becomes 


q(0 + ) 


1.15 [ .29 n (O z ) + .5 n (O)] /n (N 2 ) 

1+1.5 [n(O z ) + .5 (0)J/n (N.,) * q TOTAL 

(2.9) 


q(O^) 


1. 5 [ .71 n (O z )l /n (N 2 ) 

1 + 1. 5 In (O z ) + .5n (O)J /n (N 2 ) 


q TOTAL 

( 2 . 10 ) 


q(N 2 ) = .81 fq TOTAL - (q(O^) +q (0 2 ))] (2.11) 


q(N + ) = . 19 [ q TOTAL • (q(o + ) J q (Oj))! (2.12) 


Equations (2.9) through (2. 12) have the effect of distribute g the 
total ion- electron pair production rate q among the various ions as 
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0.62 : 0. 17 : 0. 14 : 0.07 : 1.00 for q(N*) : qfO*) : q(N + ) : q(0 + ) : q 
as reported by Swider (1969) for the E- region. 


2.3 Chemistry 

The ion neutral reactions of importance in the E- region used 
in this work were reviewed by Biondi (1969), Bortner et al. , (1972) 
and Ferguson (1974). Data was available for the ion temperature 
which was included into the electron density model. Those rates 
which varied as a function of temperature did so according to an 

_ 3 , 

inverse law (T) where a is a measured quantity. 

A reaction scheme including the more important reactions 
(Ferguson, 1974) is shown in Figure 1. The steady state of the ion 
densities is solved using the continuity equation 


d [ n .] 

”dt 


= Q(n.) - l(n.) + p(n.) - div (n.v.) 

J J J J J 


(2.13) 


where 


Q(n.) is the photoionization of species n.. 

l(n.) is the loss due to ion- neutral charge transfer. 


PU*: 


is the gain due to ion-neutral charge transfer. 


div (n. v.) is the divergence of the flux n. v.. 

J J J J 

The equation is solved for the five ions 0 + , O^, N + and NO + . While 
N* is a major ion produced because of its large concentration, it 
is not an abundant ion in this region because of its fast reactions 
with C >2 and O. Because of their low ionization potentials most 
ionization takes the form of NO + or O*. 



Figure 1 Chemical Reaction Diagram 
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2.4 Ion Transport 

Several authors (Monro, 1970; Strobel, 1974; Ogawa and 
Shimazaki, 1975) have concluded that the effects of ion drift in 
the E- region are neglegible. In this study interest was centered 
around exploring rather small departures from simple solar zenith 
angle dependencies and therefore in every minor contributing effect. 

The ion drift /elocities were measured at the Arecibo 
incoherent scatter facility in the vertical direction for the E-region 
along with the electron densities. The effect of ion drift velocities 
on the electron density was investigated by its inclusion in the 
continuity equation through the divergence term div (n^v). An idea 
of the effect of the divergence term can be obtained by assuming that 
the departures from photoequilibrium are small and that the 
dominant recombination process is dissociative and can be 
represented by a coefficient a. With this approximation at 117 km 
when the electron density vertical gradient is small the perturbation 
from a photoequilibrium density is 


An ~ 
e 


1 dv 
2 a dz 


(2.14) 


where v is the vertical velocity and at 100 km when the vertical 
gradients in velocity are small 


A n 

e 


v 1 

2 a n 

e 


d n 

e 


dz 


(2. 15) 


CHAPTER III 


DISCUSSION OF THE ANALYSIS 

3. 1 Causes of Electron Density Variations 

Electron densities measured by incoherent scatter sounding 
were compared with calculations derived from the model described 
in section 2. 1. The incoherent scatter measurements provided 
the following information: 

1. electron density, 

2. collision frequency in the region of 100 to 117 km. This 
can be used to extract the total density. 

3. temperature, 

4. vertical amplitude of the ion drift. 

An examination of model and experimental electron density 
profiles in the E-region shows very clearly two regions, a higher 
region where the electron density gradient with altitude is very small 
and a lower region of large vertical electron density gradient. 
Varying the total solar flux has the effect of multiplying the electron 
density profile by a constant factor at all altitudes. Varying the 
total density has the effect of varying the height of the profile without 
changing its shape. A multiplication factor for both the density and 
solar flux was chosen for each set of data corresponding to measure- 
ments of one day. 

This study investigated 1 the effect of variations in density, 
pressure, temperature, vertical ion drift velocity and atomic and 
molecular oxygen density on the electron density calculated from 
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Some idea of the vertical velocities in this height range can be 

obtained from Harper et al. , (1975) who found at Arecibo divergences 

- 3 - 1 

of vertical velocities at 117 km of the order of 6 x 10 sec and 
vertical velocities of less than 5 m sec * at 100 km. 

Using these formulae with the values estimated above from 
Harper et al. , (1975) wind measurements give values of A n^/n^ of 
the order of 0. 1 at 117 km and less than 0.01 at 100 km. The effect 
at 100 km is not large but at 1 17 km is not negligible compared to 
the morning to afternoon differences. 


the model described in Chapter II. The investigation was 
separated between the two regions as described previously with the 
division occurring at about 110 km. 

The sensitivity of the electron density to changes in the 
neutral density, pressure, temperature, ion drift, and relative 
atomic and molecular oxygen were investigated in the lower region. 
Sensitivity factors for density, pressure and temperature were 
determined by allowing the model to be varied by one of the three 
parameters and the resulting electron density compared with model 
results when no parameters varied. The calculations were done 
for three days, August 12, 1974, September 17, 1974 and September 
18, 1974 and the results are shown in Tables 3, 4, and 5. From 
the relative changes in the electron density due to the three parameter 
it is evident that no one parameter consistently dominates in the 
lower E-region. As discussed in section 2.4 the relative change 
in electron density at 100 km due to vertical ion drifts is on the 
order of . 0 1 . In Figure 2 the contours of electron density as a 
function of atomic and molecular oxygen are shown for 100 km. The 
electron density is seen to be almost independent of atomic oxygen. 
The maximum sensitivity of electron density to molecular oxygen 
variations is .05. 

The same analysis was performed on the upper region. Again 
not one of the three parameters, density, pressure and temperature, 
consistently dominated over the other two. As shown in section 2.4 
the vertical wind velocity accounts for a relative change in electron 
density of . 10 which is of the same order of magnitude as that due to 


Table 3: Relative Electron Density Variations for August 12, 1974 
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Altitude 

Pressure 

Temperature 

Dens ity 

Pressure 

Temperature 

Dens ity 

117 

.021 

.086 

.087 

10.4 

10.0 

12.2 

114 

.017 

.079 

. 102 

12.2 

11.6 

11.8 

110 

. 020 

. 066 

. 053 

12.4 

12.8 

12.4 

107 

. 022 

.046 

.019 

12.4 

15.2 

13.0 

104 

.051 

.087 

.037 

12.6 

15.6 

13.2 

100 

.096 

. 196 

.200 

12.6 

15.8 

13.6 


Table 4: Relative Electron Density Variations for September 17, 1974 


= A sin ( ii" - t c )) 
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A 



t 
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Altitude 

Pressure 

Temperature 

Density 

Pressure 

Temperature 

Density 

117 

.056 

.077 

. 066 

14. 0 

12. 0 

4.0 

114 

.050 

.078 

. 050 

14. 2 

12. 5 

6. 0 

110 

.033 

.089 

. 040 

6.6 

14. 0 

6.4 

107 

. 027 

. 041 

. 034 

6. 6 

17.2 

11.2 

104 

. 069 

. 006 

.098 

6. 6 

18.6 

11.4 

100 

.032 

. 170 

. 195 

11.5 

18. 0 

11.5 



Table 5: Relative Electron Density Variations for September 18, 1974 
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A ltitude 

Pressure 

Temperature 

Dens ity 

Pressure 

Temperature 

Density 

117 

. 000 

. 1 14 

.025 


12.4 

17.6 

1 14 

. 006 

. 1 16 

.033 

18.8 

13. 0 

17.8 

110 

.012 

.089 

.055 

18.6 

13.4 

19.8 

107 

.028 

.049 

.019 

19.6 

16.6 

20.6 

104 

.092 

.044 

.033 

19.6 

18.8 

18.8 

100 

. 191 

. 146 

.205 

19.6 

19.6 

19.8 
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density, pressure and temperature variations. The contours of 
electron density as a function of atomic and molecular oxygen are 
shown in Figures 2 and 3 revealing a dependence of electron density 
on both constituents. The maximum gradient or sensitivity is .05 for 
atomic oxygen and .02 for molecular oxygen. 

3.2 Feasibility of Determining Neutral Parameters from Electron 
Density Variations 

The three days, August 12, 1974, September 17, 1974 and 
September 18, 1974 were used to investigate the feasibility of 
determining variations in the neutral atmosphere from deviations 
in electron density from a basic secant X variation in electron density. 
If one of the parameters totally dominated the others in causing 
electron density variations it would be possible to work back from the 
electron density to obtain that parameter's variation. From Tables 
3, 4, and 5 it is seen that the relative changes in the electron density 
corresponding to changes in the temperature, density and pressure 
are comparable implying that the electron density variations do not 
reflect a variation in just one parameter. 

If there existed a relationship among the atmospheric param- 
eters it would be possible to determine their variations. For the 
three days examined there existed a prevalent semi-diurnal tide. 

At 100 km the phase difference between the temperature and the 
density were 2, 7 and 0 hours for August 12, 1974, September 17, 

1974 and September 18, 1974 respectively. There existed no constant 
phase relationship at any one height. 
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Figure 2 Electron Density as a Function of Atomic and 
Molecular Oxygen Densities at Noon for an 
Altitude of 100 km 
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CONTOUR LINES OF ELECTRON DENSITY (mW ) 



Figure 3 Electron Density as a Function of Atomic and 
Molecular Oxvgen Densities at Noon for an 
Altitude of 1 1 7 km 
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It would be very valuable if the amplitude and phase of the 
tidal motions could be determined from electron density measure- 
ments alone since they are available simultaneously at a large 
number of stations over a long period of time and for a large variety 
of conditions. It is the conclusion of this study that it is not possible 
to do this. 

3.3 Comparison of the Modeled and Measured Electron Densities 

The basic shape of the measured electron densities in the 
E-region varies as the secant of \ , the solar zenith angle, however 
there are deviations which can not be explained by a changing solar 
zenith angle alone. An objective of this study was to investigate the 
effect of neutral atmospheric variations on the electron density. 

As stated in section 1. 2 care must be taken in using a solar flux 
model for times other than when measured. The solar flux for one 
day for all wavelength groups was multiplied by a constant in order 
to obtain the level of ion density measured in the upper section of the 
E-region where the gradient in electron density is about zero. The 
density at all altitudes was multiplied by a constant in order to 
obtain the level of ion density measured in the upper section of the 
E-region where the gradient in electron density is about zero. The 
density at all altitudes was multiplied by a constant in order to have 
the proper altitude of the division between the upper and lower 
regions . 

The effects of measured temperature, density and vertical 
ion motion were included in the model of electron density. The 
inclusion of these variations did not consistently eliminate the 
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discrepancies between the modeled and measured electron densities. 
An attempt was made at reconciling the differences by varying the 
relative amounts of atomic and molecular oxygen. To obtain a 
good agreement a 45 per cent variation in molecular oxygen was 
required for atomic oxygen. Variations of these magnitudes are 
unexceptably large when compared with theoretical variations (Mayr 
and Volland, 1972) or measured variations (Alacayde et al. , 1972; 
Danilov, 1972 and Offerman, 1974). The difficulty in using atomic 
and molecular oxygen to reconcile discrepancies between measured 
and modeled electron densities was discussed by Monro, Nisbet 
and Stick ( 1974) . 

The basic problem in an analysis of this type is the number 
of independent variables. In the present problem the major 
uncertainties were in the relative densities of atomic and molecular 
oxygen to the molecular nitrogen density and the spectrum of the 
EUV solar flux. It would seem that further studies of this kind to be 
profitable should be done in conjunction with a series of observations 
such as those described by Torr et al. , (1976) in which satellite 
observations of these important parameters are made simultaneously 
with the incoherent scatter measurements. 

3.4 Specific Conclusions 

E- region density dependence on neutral atmospheric param- 
eters was determined for the upper and lower altitude ranges. The 
results show that variations in temperature, density and pressure 
cause variations in the electron density of comparable amounts. 



While the electron density in the lower region is sensitive to variations 
in the relative molecular oxygen density it is sensitive to variations 
in both atomic and molecular oxygen density in the upper region. The 
effect of ion drift while small in the lower region is not negligible in 
the upper region. 

The effects of tides in atmospheric parameters or the electron 
densities were studied with the objective of determining the feasibility 
of relating measured electron density variations back to one or more 
neutral atmospheric parameter variations. Since the major 
contributing parameters have comparable effects and there exists no 
fixed relationships between them, it can be concluded that tidal 
variations in the atmospheric parameters can not be determined 
from electron density variations alone. 

The modeled electron densities including the effect of tempera- 
ture, density and vertical ion drift measurements were compared with 
measured values of electron density. While some variations were 
reproduced there existed differences between them. The basic 
problem in determining the cause of discrepancies is the number of 
independent variables. It is concluded that simultaneous measure- 
ments of important parameters are necessary to sort out the causes 
of these discrepancies. 
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